Systematics of the charged-hadron Pt spectrum and the nuclear suppression factor in 
heavy-ion colhsions from ^/snn = 200 GeV to ^/s^n = 2.76 TeV 
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In this paper, our goal is to make a simultaneous analysis of the high- and low- Pt parts of the 
charged-hadron Pt spectrum measured by the ALICE collaboration [l[ in central Pb-Pb collisions 
at y/SNN = 2.76 TeV at the Large Hadron Collider (LHC), based on models which have been 
successfully applied and constrained in Au-Au collisions at Relativistic Heavy Ion Collider (RHIC). 
For the hydrodynamical modeling with which we obtain the 1ow-Pt spectrum, we have computed 
the initial conditions based on perturbative QCD (pQCD) minijet production and saturation. The 
sensitivity of the obtained charged-hadron Pt spectrum on the hydrodynamic model parameters 
is studied. For the high-Pr part, we apply a number of parton-medium interaction models, which 
are tuned to describe the nuclear suppression factor Raa measured at the RHIC in central Au-Au 
collisions at y^Siviv = 200 GeV. We find that the higher kinematic reach of the LHC, manifest in 
the hardening of the pQCD parton spectral slope, is in principle very efficient in discriminating the 
various models. However, due to the uncertainties in the p-p baseline, none of the tested models can 
be firmly ruled out with the present ALICE data. Comparison with the LHC data in this approach 
shows that the combined hydrodynamic and pQCD4-jet quenching components can reproduce the 
data well in the whole measured Pt range. 

PACS numbers: 25.75.-q,25.75.Gz 
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I. INTRODUCTION 



The very first ALICE physics results from the ^^snn = 
2.76 TeV Pb-Pb collisions at the LHC were the ellip- 
tic flow ^1, multiplicities Q, as well as the transverse 
momentum (Pt) spectrum of charged particles In 
the measured charged-particle Pt spectrum, which cov- 
ers the range Pt < 20 GeV and extends over three orders 
of magnitude, two quite distinct regions can be seen for 
the centralmost collisions: at Pt ^4 — 5 GeV the spec- 
trum is exponential, while at Pt ^ 5 GeV it shows a 
power-law like behaviour. Both of these carry very im- 
portant information on the QCD dynamics of the sys- 
tem: the \ow-Pt spectrum, responsible for the hadronic 
bulk multiplicity, reflects the collective transverse motion 
(flow) developed in the system during its entire spacetime 
evolution, while the high-Py spectrum, and its suppres- 
sion relative to the yield in proton-proton collisions, tells 
us about energy losses of high transverse momentum (pt) 
partons on their way out of the dense bulk matter. 

In this paper, our goal is to make a simultaneous anal- 
ysis of the high- and \ow-Pt parts of the charged-hadron 
Pt spectrum measured by the ALICE collaboration [1,] 
in Pb-Pb collisions, based on models which have been 
successfully applied and constrained in Au-Au collisions 
at RHIC. While the extrapolation of the pQCD parts 
from y/SNN = 200 GeV to 2.76 TeV is straightforward. 
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the nonperturbative features in bulk dynamics and its 
connection to partonic energy loss contain considerable 
uncertainties. 

A generally accepted framework for describing the bulk 
QCD-matter evolution in A-A collisions, is relativistic 
hydrodynamics, which at the same time is the only frame- 
work which can incorporate the effects of the phase tran- 
sition predicted by lattice QCD. There has been consid- 
erable progress in the hydrodynamical modeling of ultra- 
relativistic heavy-ion collisions over the last years, and 
into many directions, too: one has moved from solving 
1-1-2 dimensional boost symmetric ideal hydrodynamical 
equations to genuinely 1-1-3 D ideal hydrodynamics 
[IHI, and to l-h2 D and even 1-^3 D dissipative 

hydrodynamics p^ . Furthermore, various hybrid models 
have recently been developed, where the hydrodynamical 
evolution is coupled with a hadronic cascade afterburner 
[1, [T^ - [l7| . Also genuine event-by-event hydrodynamic 
studies have been performed, both in the ideal hydrody- 
namics 11814231 as well as in the viscous case uM- 
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Hydrodynamics itself does not provide any prescription 
for the extrapolation in y^SNN- The ^/snn dependence 
of bulk dynamics enters mainly via the initial conditions 
(IC) which are crucial element in these studies. One can 
try to fit the IC using the available data (multiplicities, 
Pt spectra, elliptic flow) as constraints. Then, however, 
the more fltting is done the more one looses in predictive 
power. To improve upon this, the IC may be computable 
using a dynamical model for primary production of QCD 
quanta, which provides the needed ^/snn dependence. 
One such initial-state model, which we will employ in this 
paper, is the pQCD -I- saturation approach, the "EKRT 
model" [13 . When combined with ideal hydrodynam- 
ics, this model correctly predicted the charged-hadron 
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multiplicities in y/sNN — 56, 130 and 200 GeV Au-Au 
collisions at RHIC, and within 7 % also at the LHC [23| . 
Also the centrality dependence of the multiplicity is con- 
sistent with the data, see Fig. 23(a) in [26] and Fig. 4 
in [l^l , and the measured Pt spectra of the bulk of the 
identified hadrons at RHIC have been reproduced well in 
this framework [2a|. Elliptic flow has also been success- 
fully addressed [H,!!^, and in particular, the similarity 
of the differential elliptic flow at RHIC and LHC as well 
as the increase of integrated elliptic flow from RHIC to 
LHC was predicted in [28(|. Furthermore, the emergence 
of the pQCD tail from the hydrodynamic spectrum at 
Pt ~ 4... 5 GeV in central Pb-Pb colhsions at the LHC 
was predicted in [25| . 

High-pT partons are created in hard pQCD subpro- 
cesses along with bulk multiplicity production in ultra- 
relativistic heavy-ion (A-A) collisions. The idea that the 
final state interaction of such partons with the surround- 
ing matter reflects properties of the bulk and that hence a 
measurement of high Pt hadrons can be used as a probe 
of the QCD medium, is known as "jet tomography" [sol - 
|35| . One of the expected signatures of this final state 
interaction is the suppression of high Pt hadron produc- 
tion in A-A collisions as compared to the scaled expec- 
tation from proton-proton (p-p) collisions, measured by 
the nuclear suppression factor Raa ■ This phenomenon is 
often called "jet quenching" (although this is is slightly 
misleading since the observable is not a jet of hadrons 
but inclusive single-hadron spectrum). Experimentally, 
jet quenching has been measured at RHIC with great 
precision as a function of collision centrality and orienta- 
tion with respect to the reaction plane [36| . 

A number of parton-medium interaction models has 
been proposed and tested against the experimental re- 
sults for Raa together with a well-constrained hydro- 
dynamical description of the bulk medium [37H39| . How- 
ever in systematical comparisons of models with the data, 
even including observables such as high Pt back to back 
correlations, ambiguities remain |40l442l | which do not 
allow firm answers as to what the correct dynamics of 
parton-medium interaction is in nature. A partial answer 
was obtained with regard to models proposing elastic (or 
more general incoherent) energy loss of leading partons 
[431 - 1471 . Such models fail to reproduce pathlength de- 
pendent observables such as the spread between in-plane 
and out of plane emission of hard hadrons [3 113 ■ This 
result establishes that quantum coherence is a crucial in- 
gredient of the dynamics, but fails to discriminate among 
the various QCD radiative energy loss models to the de- 
gree that a quantitative extraction of medium parameters 
would be possible with good accuracy [50| . 

The underlying reason for this failure is that the nu- 
clear suppression factor is largely independent of the 
functional form of the energy loss probability distribu- 
tion P{AE) for a given leading parton 51]. While differ- 
ent models predict different forms for P{AE), the pQCD 
parton spectrum at RHIC kinematics is so steeply falling 
that even a small shift in parton energy due to energy loss 



to the medium effectively acts like a complete suppres- 
sion of the parton. Thus, only a small part of P{AE) 
close to zero energy loss is actually probed in observ- 
ables. The vastly larger kinematic range of the LHC, 
leading to a significantly harder pQCD parton spectrum 
is expected to change this situation and to allow to probe 
more deeply into P{AE). 



While leading-parton energy loss models are sufficient 
to describe single inclusive high Pt hadron production, 
there is a second class of models which treat the whole 
medium-modification of a parton shower in the medium 
,52-56] with the aim of eventually describing fully recon- 
structed jets in heavy-ion collisions. These models are 
often Monte-Carlo (MC) codes extending vacuum shower 
codes such as PYTHIA ^] or HERWIG ^] to include 
the interaction with a medium. A systematic comparison 
of these codes with pathlength dependent observables is 
so far largely absent. 



In this paper, we make a simultaneous analysis of 
the high- and low-Py parts of the charged-hadron 
Pr spectrum measured by the ALICE collaboration 
[i[ in ,/sNN ~ 2.76 TeV central Pb-Pb collisions 
at the LHC. We first study the compatibility of the 
pQCD-|-saturation-|-hydrodynamics framework with the 
measured low-P^ spectrum, systematically charting the 
model uncertainties. Then for the high-P^ part, using 
the obtained hydrodynamical evolution of the bulk mat- 
ter as background, we investigate what discriminating 
power the first measurement of Raa at the LHC [l| of- 
fers for models which are tuned to describe the observed 
nuclear suppression at RHIC. We consider both leading- 
parton radiation 159|1 and elastic energy losses as 
well as showers [53, [SJ] modified by the QCD medium. 
For other recent works discussing Raa at the LHC, see 
Refs. [63,[6l- 



For the Raa study, our strategy is as follows: First, 
in moving from RHIC to LHC using our default hydro- 
dynamical set-up, we compute Raa for different models 
of parton-medium interaction without any re-tuning of 
parameters (straight extrapolation). Since some amount 
of uncertainty is expected to originate from the uncer- 
tainties in the hydrodynamical initial state, and due to 
the fact that energy-loss model parameters are known 
to differ for the same parton-medium interaction model 
even among constrained hydrodynamical models (42j , we 
re-tune in a second run the model parameters to the best 
fit of the ^fswN = 2.76 TeV data and quote the difference 
to the tune at ^/snn = 200 GeV. This allows to gauge 
how the change in ^/snn acts to constrain models. Using 
the best fits which we obtain for the 1ow-Pt and high- 
Pt parts, we investigate to what extent the measured 
charged-hadron Pt spectrum can be reproduced. 
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II. THE HYDRODYNAMICAL BULK 
DESCRIPTION 

A. pQCD+saturation+hydrodynamics framework 

For obtaining the bulk hadron px spectra in a hydro- 
dynamical framework, we need to define the QCD mat- 
ter initial conditions, solve the hydrodynamical equations 
numerically with a given Equation of State (EoS), com- 
pute the thermal particle spectra at freeze-out, and ac- 
count for the strong and electromagnetic decays of un- 
stable particles. 

We compute the QCD matter initial densities using 
the EKRT saturation model 23], which is based on 
collinearly factorized pQCD minijet production and the 
conjecture of gluon saturation in the transverse plane. 
In this model, which has been quite successful in pre- 
dicting the multiplicities from RHIC to LHC, saturation 
of primary parton (gluon) production is assumed to take 
place when the minijet production vertices, each of which 
occupies a geometric uncertainty area ^ t^/Pt^ start to 
overlap. For central A-A collisions studied here, the fol- 
lowing criterion is fulfilled at saturation: 



(1) 



where the number of produced minijets at pr > po in the 
mid-rapidity unit Ay = 1 can be written in terms of the 
standard nuclear overlap function Taa and leading-order 
(LO) pQCD cross sections as N(j)q) = TAAiO)a{N), with 
(see Ref. ^23] for details) 
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Above, TAA{h) = / d^sTA{s)TAis - b), where TA(r) = 
J dzuAir) is the nuclear thickness function computed 
from the Woods-Saxon nuclear densities ua with ng = 
0.17 fm^^, d = 0.54 fm. The inclusive cross section for 
producing partons of flavours k and I above is, as usual. 



^^AA^kl+X 

dpldy,dy2 =KY.^ih/A{x,,Q')x2l„A{x2.Q') 



da'3- 



dt 
(3) 

For the nuclear parton distribution functions (nPDFs), 
we use the CTEQ6L1 PDFs ^] together with the EPS09 
nuclear effects [63|. For transverse energy (Et) produc- 
tion, the pQCD minijet calculation can be extended to 
next-to-leading order (NLO) in an infra-red and collinear 
singularity safe manner (whereas the number of produced 
minijets is not well defined beyond LO) |6J|]. The corre- 
sponding updated K factors are, however, not yet avail- 
able [65| . which is why we simply perform a LO pQCD 
calculation here and fit the K parameter in Eq. ([3]) so 
that the minijet production at saturation leads to (af- 
ter hydrodynamic evolution and resonance decays) the 



measured charged hadron multiplicity in the LHC Pb- 
Pb collisions at y^SNN = 2.76 TeV in the 0-5% centrality 
class, dNch/dri = 1584 [3]. The obtained K thus also ac- 
counts for all higher-order contributions. The centrality 
selection is simulated by considering a central ^off-^off 
collision of an effective nucleus ^eff = 192 < A, as ex- 
plained in Ref. [53] . 

Most importantly, with the EKRT-based modeling 
above, we can also estimate the formation time of the 
system from the pQCD dynamics: according to the un- 
certainty principle, t; 1/psat, where psat is the solution 
of Eq. III). With K = 1.54, obtained by an iterative fit 
to the measured multiplicity, we have psat = 1-58 GeV 
and Ti = 0.12 fm. 

In the EKRT framework, where the computation of 
the minijet Et production is possible in NLO, the initial 
conditions should be given in terms of the energy density 
instead of the entropy density. Once the saturation scale 
is known, we compute the local initial energy density by 
distributing the minijet Et over the transverse plane at 
the time according to (for more details, see |25| ) 



e(r,Tj) = TA{r)TAir) 



<j{Et) 
nAy 



where 



a{ET) = / dpTPT 



dp'^ 



(4) 



(5) 



Above, we distribute the energy density into the trans- 
verse plane according to the binary collision profile, thus 
our default set-up is the "eBC" initialization. It should 
be emphasized, however, that since we do not attempt 
to make the saturation condition ([ij strictly local in the 
transverse plane (which would lead to a varying psat , see 
(27I [29| for such discussion) , the transverse profile is not 
uniquely fixed here. 

We use ideal hydrodynamics to describe the space-time 
evolution of the bulk matter. Since we will consider only 
mid-rapidity observables here, longitudinal boost invari- 
ance as well as neglecting the net-baryon number are very 
valid approximations. We solve the 2-1-1 dimensional rel- 
ativistic hydrodynamical equations, dfj.T^^" = 0, using 
the SHASTA algorithm J^^. For the EoS which closes 
the set of dynamical equations, we choose the recently de- 
veloped EoS s95p-vl from Ref. [gl]. We assume a very 
rapid thermalization here, taking the formation time ti 
as the starting time tq for the hydrodynamical evolution. 
Furthermore, a full chemical equilibrium and zero initial 
transverse flow are assumed. 

Regarding these initial conditions, we should empha- 
size three points here. First, the early initialization of 
the hydrodynamical evolution, tq oc l/jOsati is quite es- 
sential, since, as pointed out long ago [60], pQCD mini- 
jets do produce a large amount of Et, and only by doing 
PdV work over a long enough time early on, the hydro- 
dynamically evolving system can sufficiently degrade its 
transverse energy: As shown in [241 ] . the amount of the 
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measured final-state Et is only a third of the initially pro- 
duced Et- Second, although the system may not be fully 
thermal at early times, the early start accounts for the 
buildup of flow and pressure as well as PdV work during 
the thermahzation stage. Third, as discussed e.g. in [70| . 
thermal photon production is very sensitive to the hydro 
initialization time. In order to explain the photon pro- 
duction measured at pr ~ 2 GeV in ^/sNN — 200 GeV 
Au-Au collisions, one needs a substantial thermal pho- 
ton production component from the QGP. This can be 
obtained only if the initialization time is small enough. 
For these reasons, we believe that to start the hydrody- 
namical evolution at an early time with zero transverse 
flow is physically a well motivated approximation. 

Finally, hadron spectra in the hydrodynamical model 
are calculated from a constant temperature freeze-out 
hypersurface using the Cooper-Frye formula [71] , and ac- 
counting for the strong and electromagnetic decays of un- 
stable particles. The freeze-out temperature Tp is flxed 
so that we can describe the measured positive pion spec- 
tra in ^/snn = 200 GeV Au-Au coUisions at RHIC. The 
value of Tp is found to be 165 MeV for the eBC pro- 
flle. As shown in [72], for computing the Pt spectra of 
hadrons, keeping the Tp unchanged from RHIC to LHC 
is a good approximation to a more dynamical decoupling 
treatment where the scattering rates are compared with 
the expansion rate of the system. 



B. The results: Hydrodynamical pr spectra and 
their systematics 

Figure [T] shows the hydrodynamically obtained px 
spectrum of charged hadrons in 0-5% central ^s^vtv — 
2.76 TeV Pb-Pb collisions at the LHC, and its compar- 
ison with the ALICE data [ll]. Also shown is the com- 
parison with the computed pQCD-part of the spectrum, 
which is subjected to the energy losses of high-energy par- 
tons discussed in Sec. UTT] below. As seen in the flgure, 
the agreement of the hydrodynamically obtained spec- 
trum is quite good down to 4-5 GeV, where the pQCD 
tail takes over. The observed behaviour, the change from 
the hydro-dominated to the pQCD-dominated spectrum 
at = 4 — 5 GeV is indeed very similar to what was 
predicted in 25| (see Fig. 15 there). 

To study the sensitivity of the hydrodynamic pr spec- 
trum to the model uncertainties, we perform the follow- 
ing systematics shown in Fig. [21 

First, in the upper left panel, we vary the freeze-out 
temperature between 120 MeV and 165 MeV, keeping 
the initial conditions fixed to the default set-up. With 
the single- scenario we use here (i.e. no partial chemi- 
cal equilibrium or coupling to a hadronic afterburner), a 
lower Tp enables more transverse flow to develop which 
in turn leads to a flatter Pt spectrum which easily over- 
shoots the data both at RHIC and at the LHC. Thus, 
with the eBC and small- tq set-up we need a high value 
of Tp. 
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FIG. 1: (Color online) The transverse momentum Pt spec- 
trum of charged hadrons in 0-5% central ^/snn ~ 2.76 TeV 
Pb-Pb collisions as measured by the ALICE collaboration 
[3| and compared with theoretical calculations using a two 
component picture: the low Pt region is described by 
pQCD-f saturation+hydrodynamics whereas in the high Pt 
region we apply a pQCD + jet quenching picture, here in the 
ASW formalism (see text). Shown for comparison is also the 
pQCD result without jet quenching. Due to the Pt depen- 
dence of the jet quenching, this has a different shape which 
would not agree with the data. 



Second, the upper right panel shows the sensitivity of 
the results to the initial time: here tq has been var- 
ied in the range 0.5/psat ~ 2/psatj keeping the initial 
minijet Et and the final Tp unchanged. We observe 
that the eBC set up favours a small hydro initializa- 
tion time but that at the same time - since the ini- 
tial Et is conserved and not the initial entropy (^o) 
- the multiplicity decreases slightly with decreasing tq: 



ch 



So - e^/^To ^ (i?T/To)=^/Vo 



e: 



3/4 1/4 



Thus, 



changing tq by a factor 2 causes a change of 19 % in the 
multiplicity - a change which is already well beyond the 
5 % error bar in the data 3 . 

Third, the panel on lower right shows the sensitivity 
of the hadronic Pt spectrum to the choice of the energy 
density transverse profile. In addition to the binary pro- 
file TA{r)TA{r) in Eq. ([1]), we have distributed the minijet 
Et also over a wounded nucleon density profile (eWN), 
keeping however tq fixed. In this change, the charged 
particle multiplicity then increases from 1580 to 1640. 
Based on fitting the measured charged-particle Pt spec- 
trum in y^SNN = 200 GeV Au-Au collisions at RHIC, we 
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FIG. 2: (Color online) Systematics of the hydrodynamic Pt spectrum in 0-5% central ^/sNN = 2.76 TeV Pb-Pb collisions. 
Keeping the initial energy fixed, we show the sensitivity of the computed spectrum to the freeze-out temperature Tp (upper 
left panel), to the hydrodynamic initialization time to (upper right), and to the initial transverse profile (lower right). Allowing 
for a change of 20% for the computed multiplicity, the sensitivity of the results to the value of K (affecting both the initial Et 
and to) is shown in the lower left panel. The data points are from the ALICE measurement [l|]. 



use Tp = 160 MeV for the eWN case. Due to the slower 
build-up of transverse flow in the eWN case, the eWN 
spectrum becomes slightly steeper than our default eBC 
case with Tp — 165 MeV. From this, we can deduce that 
fine-tuning the transverse profile towards an eWN profile 
would be certainly possible but also that an early initial 
time To ~ 1/psat is still required. 

Fourth, the lower left panel shows the sensitivity of the 
spectrum to the fit parameter K. For this figure, we al- 
low a change of 20 % (ca. 300) in the multiplicity, and 
compute the Et and tq separately in each case but keep 
Tp constant. We can see how an increase(decrease) in the 
multiplicity corresponds to a larger increase(decrease) in 
the value of K but a smaller decrease(increase) in tq. 
These scalings can be deduced from Eqs. (l)-(5), keep- 
ing track of the powers of po and K: At the scaling limit 
a{po) ~ K/pI, and Eq. ^ gives ftat - K^'^- The en- 
ergy density scales as e o-(£^T)ftat ^ (^/Psat)Psat ^ K, 
and the charged-particle multiplicity then scales as iVch ^ 
So ^ f?I^To ^ K^l'^. Thus, a 20% increase in the multi- 
plicity corresponds to a 40 % increase in the value of K 



and a 10 % decrease in tq. We see on the one hand, that 
our initial-state modeling is fairly robust against small 
variations in K and on the other hand that fine-tuning 
of K would depend on the transverse profile and also on 
Tp- Such fine-tuning, however, we do not consider very 
meaningful before more data, on e.g. identified hadron 
spectra, are available. 



From Fig. [21 we conclude on the one hand that with our 
pQCD-|-saturation-|-(ideal)hydrodynamics framework we 
are committed to a fairly narrow window of the param- 
eters Tq and Tp, and on the other hand that the profile 
uncertainty can be considered to be the main uncertainty 
from the jet quenching viewpoint. Further fine-tuning of 
the hydrodynamical description is left as future work, 
when we are studying the centrality dependence of the 
hadronic Pt spectrum. 
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III. PARTON-MEDIUM INTERACTION 
MODELS 

The starting point for a computation of the high Px 
hadron yield in an A-A collision is the initial spectrum 
of hard partons. In our framework, the differential cross 
section da^^^^^^-^ for the production of a parton / in 
an A-A collision is calculated in LO pQCD by integrating 
out the unobserved parton kinematics in Eq. ([3]) (explicit 
expressions are given in [40| and references therein). 

Uncertainty relation arguments indicate that the 
medium cannot modify the hard process itself, but rather 
influences the fragmentation pattern of the outgoing 
highly virtual partons, in particular their development 
into a parton shower. Thus, for in-medium shower 
models, the expression to evaluate is the convolution of 
the partonic production cross section with the medium- 
modified fragmentation function (MMFF), 



computed from the convolution 



""mod 



AA^f+X 
vac 



(6) 



where {D^^^j^{z, ^'^))taa is the MMFF, averaged over the 
geometry of the medium, z is the fractional momentum 
of produced hadrons given a parton /, and /i^ is the 
hadronic momentum scale. To evaluate this expression 
requires knowledge of both the geometry of the medium 
(e.g. in terms of a spacetime description of medium den- 
sity) and of the MMFF Dmm{z, fJ-^, C) for any given path 
^ through the medium. 

If the angle between outgoing parton and the reac- 
tion plane is </>, the path of a given parton through the 
medium C('''), i.e. its trajectory ^ as a function of proper 
medium evolution time r is determined in an eikonal ap- 
proximation by its initial position ro and the angle (j) as 
C(t) = (xq + T cos((/)), 2/0 + Tsin((/))) where the parton is 
assumed to move with the speed of light c = 1 and the 
x-direction is chosen to be in the reaction plane. How 
Dmm{z, ^pX) is obtained once a medium is specified is 
characteristic for a given model of parton-medium inter- 
action and will be discussed for the code YaJEM later. 

Once the MMFF for a given path is known, the aver- 
aging over the medium geometry is given by 



{Dmm{z,IJ^p))taj 



^ l'2'K nOO nOQ 

— dcj) dxo dyoPixo,yo)DMMiz, iJ-l, 

-^TT J -oo J -oo 



C). 



(7) 



Here, the initial distribution of hard vertices in the trans- 
verse {x, y) plane is assumed to be calculable as 

p, ^ T^(ro+b/2)T^(ro-b/2) 

In the leading-parton energy loss approximation, the 
medium-modified production of high-P^ hadrons can be 



da. 



AA^h+X_ 
med 



E 



dati^f+^^{P{AE))TAA^Df^'^{z,^^') 

(9) 

where {P{AE))taa is the medium-induced energy loss 
probability, averaged over the medium geometry and 
Df^^'{z,fi^) is the vacuum fragmentation function for 
the production of a hadron h from a parton /, at 
fractional momentum z and hadronic momentum scale 
/i^. The underlying assumption is that the dynam- 
ics of parton-medium interactions can largely be cast 
in terms of a shift in leading parton energy and that 
hence the MMFF can be approximated by the convolu- 
tion of an energy loss probability with the vacuum frag- 
mentation function, {d{jI'}{z, h'^))taa = {P{^E))taa ® 
Df^^{z, fi"^). In this case, the energy-loss probability 
for a given path C of a parton through the medium, 
P{/S.E, C), is the ingredient to be computed within a spe- 
cific model of parton-medium interaction. 

If P{AE, C) is known, the geometrical averaging in- 
volves as above integrating over all possible initial ver- 
tices {xo,yo) with the weight of P{xo,yo) and over all 
possible orientations (f) as 



2 /'2'n" /.oo /.oo 

{P{AE))TAA^7r / / / dyoP{xo,yo)P{AE, 

(10) 

In all cases, the nuclear modification factor is com- 
puted with the given medium-modified yield of hard 
hadron production as 



C). 



RAA{PT,y) = 



dN^A/dpTdy 
TAA{h)daPP/dpTdy' 



(11) 



The details of the parton-medium interaction model 
are thus in either the energy loss probability distribu- 
tion P{AE, () for leading parton energy loss models or 
the MMFF Dmm{z, Hp, C) for in-medium shower models, 
given a specific path through the medium. In the follow- 
ing, we formulate three different types of models which 
we apply to the ALICE data. 



A. Armesto-Salgado- Wiedemann (ASW) formalism 

The detailed calculation of P{AE, Q follows the Baier- 
Dokshitzer-Mueller-Peigne-Schiff (BDMPS) formalism 
for radiative energy loss [31] using quenching weights as 
introduced by Salgado and Wiedemann [soj . 

In this framework, the energy loss probability 
P{AE, C) for a path can be obtained by evaluating the 
line integrals along ({t) as 



Lo,{ro,<P)= dCCm) and (gi) (ro, ^) = / ^(^(C) 
Jo Jo 

(12) 

with the relation 

q{() = Knicd ■ 2 • e^/''(^)(coshp — sinhpcosa) (13) 
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assumed between the local transport coefficient q{(^) 
(specifying the quenching power of the medium) , the en- 
ergy density e and the local flow rapidity p with angle 
a between flow and parton trajectory (tJ, [ill- ^mcd 
is the adjustable parameter in this framework. It is 
naturally expected to be 0(1), but in fits to data at 
y/SNN = 200 GeV the parameter takes (dependent on 
the precise hydrodynamical model) values ranging be- 
tween 3 and 10 (the latter number occurs for viscous 
hydrodynamics where the initial entropy density is lower 
than in the ideal case, see jl^l)- 

Using the numerical results of [501 and the defini- 
tions above, the energy loss probability distribution 
given a parton trajectory can now be obtained as a 
function of the initial vertex and direction {ro,(f)) as 
P(Ai;;We(r,0),i?(r,0)) = PiAE,C) for Wc(C) and R = 
2ujc{Cfl{qL{C)). In practical terms, {P{AE))ta^ is 
characterized by a fairly large discrete escape probability 
without energy loss and a very broad distribution of en- 
ergy loss ranging up to O(IOO) GeV at RHIC conditions 
(for explicit figures, see e.g. [4Ci|). 



B. YaJEM (Yet another Jet Energy-loss Model) 

The MC code YaJEM is based on the PYSHOW code 
[js'l which is part of PYTHIA fE^. It simulates the evo- 
lution from a highly virtual initial parton to a shower of 
partons at lower virtuality in the presence of a medium. 
A detailed description of the model can be found in 
[53, 54]. 

The parton shower developing from a highly virtual 
initial hard parton in this model is described as a series 
of 1 — >■ 2 splittings a — 6c where the virtuality scale de- 
creases in each splitting, i.e. Qa > Qb, Qc and the energy 
is shared among the daughter partons 6, c as = zEa 
and £'c = (1 — z)Ea- The splitting probabilities for a 
parton a in terms of Qa, Ea are calculable in pQCD and 
the resulting shower is computed event by event in a MC 
framework. In the presence of a medium, the main as- 
sumption of YaJEM is that the parton kinematics or the 
splitting probability is modified. In the RAD (radiative 
energy loss) scenario, the relevant modification is a vir- 
tuality gain 



(14) 



through the interaction with the medium during the par- 
ton lifetime Ta- This modification leads to an increase in 
radiation. In order to evaluate Eq. (|14l) during the shower 
evolution, the momentum space variables of the shower 
evolution equations need to be linked with a spacetime 
position in the medium. This is done via the uncertainty 
relation for the average formation time as 



I y _ Eb Eh 



and randomized splitting by splitting by sampling Tf, from 
the distribution 



P{Tb) = exp 



n 

in) 



(16) 



The evolution for any given parton in the shower evo- 
lution is terminated as soon as the parton reaches a min- 
imum virtuality scale Qq. The result of the partonic evo- 
lution in terms of a shower of low virtuality partons is 
then passed on to the Lund model [tI] to hadronize. The 
fractional longitudinal momentum distribution of the re- 
sulting hadron distribution corresponds to the MMFF of 
the various hadron species. 

In the default version of YaJEM, the minimum virtual- 
ity scale is fixed at Qo = 0.7 GeV. In the version YaJEM- 
D (dynamical computation of Qo) [77j . the formation 
length of the in-medium shower is forced to be within 
the medium length. This corresponds to the choice 



(17) 



which depends on both in-medium pathlength L and 
shower-initiating parton energy E. The original motiva- 
tion for this prescription was to introduce a pathlength 
dependence that can account for the experimentally ob- 
served split between in-plane and out of plane emission 
of high Pt hadrons in non-central collisions [SQ] . How- 
ever, together with the stronger pathlength dependence, 
YaJEM-D also predicts a strong rise of Raa with Pt 
in angular averaged observables which we aim to test 
against the ALICE data. 

In principle, the full functional form of q{(^) could de- 
termine the MMFF, which would be computationally 
very expensive as a full MC simulation would be needed 
for every possible path in the medium. However, due to 
an approximate scaling law identified in [53], it is suffi- 
cient to compute the line integral 



AO?ot 



(18) 



(15) 



in the medium to obtain the full MMFF Dmm{z, li^X) 
from a YaJEM simulation for a given eikonal path of 
the shower-initiating parton, where fjl^ is the momentum 
scale of the shower initializing parton. The scaling law 
implies that the MC simulation has to be run only for a 
finite set of paths and makes a numerical solution of the 
geometry averaging possible. 

Each YaJEM run determines the MMFF for a fixed 
partonic scale fXp. To account for the scale evolution 
of the MMFF, runs for different fip need to be done. 
For technical reasons having to do with numerical perfor- 
mance, we like to evolve the MMFF for given hadronic 
scale as indicated in Eq. ©. For matching a partonic 
scale for which the MMFF is computed to a hadronic 
scale, we use the following procedure: For each avail- 
able partonic scale, {Dmm{z, ij,p))taa is computed, and 
the exponent n of a power law fit to the parton spec- 
trum at scale /Zp is determined. The maximum of 
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z^~^{Dmm{z, fip))TAA corresponds to the most likely 
value z in the fragmentation process, and thus the par- 
tonic scale choice is best for a corresponding hadronic 
scale Pt = z^p. The Pt dependence of the hadronic Raa 
is then computed by interpolation between runs with dif- 
ferent scale choice for the MMFF. 

As in the previous case, i^mod in Eq.(|13p serves as the 
adjustable parameter of the model once Qq is chosen. 
YaJEM requires, dependent on the underlying hydrody- 
namical model, rather natural values for i^mcd ranging 
from 0.6 to 2. 



C. Parametrized elastic energy loss 

In [S], a phenomenological model for elastic energy 
loss, consisting of a discrete parton escape probability 
and a Gaussian parametrization for the energy loss prob- 
ability was introduced to explore the pathlength depen- 
dence of incoherent energy loss. While the model it- 
self is rather simplistic, its main findings with regard 
to pathlength dependence have been confirmed later in 
a detailed MC simulation of elastic energy loss [i^. It 
needs to be stressed that unlike the previous models, the 
parametrized elastic energy loss is not meant as a seri- 
ous QCD based model for the underlying dynamics of 
parton-medium interaction. The reason that it is pre- 
sented here is rather that the simple and adjustable form 
of P{AE, C) allows insight into how the observed rise of 
Raa with Pt depends on the underlying functional form 
of the energy loss probability density, which is much less 
transparent in the context of ASW or YaJEM. 

In the model, the escape probability of a parton i with- 
out any medium interaction is computed as 



the fluctuations around the mean will grow like y/N. As- 
suming a Gaussian distribution, this leads to the ansatz 



exp 



-const. ■ a 



el 



PAiiOdC 



exp[-7i • k\ 



(19) 

where it is assumed that aei is approximately indepen- 
dent of C, and K is defined as 

K = I d^e^/*{^){coshp{0 ^smhp{^)cosa) (20) 



taking into account the flow corrections to the probed 
density. Here is a parameter with dimensions of a cross 
section measuring the interaction strength, and hence 
7g = 9/47q must hold to account for the different color 
factors of quarks and gluons. 

If the parton does not escape without energy loss, it 
must undergo a shift in energy (there is also the possi- 
bility that a strong shift into a thermal regime occurs, 
which is equivalent to an absorption of the parton) . It is 
assumed that the mean value of the shift in energy will 
grow linearly in the number of scatterings N as 

dAE = AEiaiiPMd^ (21) 
with AEi the mean energy loss per scattering, whereas 



P\AE,O^PLJ{^E)+M,exp 



{AE - a^nf 



(22) 



where Mi is a normalization such that 
d{AE)P''{AE,C) = 1 and (1121) has to hold for 
quarks and gluons separately due to the different color 
factor, ai is a parameter with the dimensions of a cross 
section times the energy shift per reaction. 

The model is thus characterized by three parameters: 

• ai controls the mean shift in energy per expected 
scattering 

• /3i governs the strength of fiuctuations around this 
mean shift. If /3i is small, the model will have a 
strong correlation between path (and hence initial 
vertex) and shift in energy, if the parameter is large, 
this correlation is lessened 

• 7i finally determines the magnitude of the escape 
probability. 

In [4^ it was discussed that the space of all possi- 
ble (Q!i,/3i,7i) which can describe the measured Raa for 
y^SNN — 200 GeV central Au-Au collisions is triangular 
and ordered by 7 — if Pesc is of the order of the measured 
Raa already, the space of allowed shifts in energy is very 
constrained. On the other hand, if Pcsc is small, more 
possibilities for a shift arise. 

Here, we investigate two distinct scenarios, one with 
large escape probabilities for quarks and gluons close 
to the allowed limit where Pg^ = 0.218, PJ',^ = 0.054 
and one with about half these values where P^^^ = 
0.12, Pg^ = 0.027. The parameters are found in Table HI 
Note that the parameters here do not correspond to the 
sets given in [48| since we use a different hydrodynamical 
model to do the extrapolation from ^/snn — 200 GeV to 
2.76 TeV here. 





small Pcsc 


large Pcsc 




0.5 


2.0 




20.0 


80.0 


7. [GeV"^] 


0.4 


0.28 



TABLE I: Parameters for the two different elastic energy loss 
scenarios described in the text. 



IV. RESULTS 

A. Pt spectra 

In Fig.[TJ we show in addition to the good description of 
low Pt bulk matter by hydrodynamics the quality of the 
description of the high Pt part above 6 GeV by pQCD 
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+ jet quenching. We conclude from this comparison that 
a LO pQCD calculation, supplemented by a JsT factor, 
is sufficient for the accuracy required to compute Raa 
reliably enough. The residual small deviations are not a 
crucial issue, since due to the nature of Raa as a ratio 
of Pt spectra deviations in the spectral shape between 
computation and data cancel to first order and only lead 
to subleading corrections. 

We also note at this point that the effect of jet quench- 
ing cannot be cast into the form of a constant downward 
shift of the unmodified spectrum, but that rather its Pt 
dependence is crucial to describe the spectral shape cor- 
rectly. As an interesting side remark, the fact that pQCD 
and jet quenching offer a good description of the spec- 
trum from 6 to 20 GeV suggests that there is fairly lit- 
tle room for another Py-dependent component of hadron 
production, such as suggested by, e.g. certain recombi- 
nation models (see e.g. |78|) in this region. On the other 
hand, e.g. the sudden recombination model ex- 
pects to see effects largely at lower momenta. However, 
since recombination models generically expect quit e dif- 
ferent effects for mesons and baryons [tI, [tI, |8]| , the 
charged hadron spectrum is not a suitable observable to 
gauge the importance of recombination in this momen- 
tum window and a detailed discussion should be based 
on identified hadron spectra. 



B. Direct extrapolation 

In Fig. [3] we show Raa computed in the various models 
for parton-medium interaction discussed above in com- 
parison with the ALICE data f]\, with their parameters 
adjusted to 0-10% central ^/snn — 200 GeV Au-Au col- 
lisions at RHIC. The assumption underlying this extrap- 
olation is that the hydrodynamical model for the bulk 
matter can be extended from RHIC to LHC in a well 
controlled manner. This is a non-trivial issue, as it is 
known that changing the hydrodynamical model at the 
same energy may amount to 50% change in model param- 
eters if the underlying dynamics is sufficiently different, 
even if both models are constrained by bulk observables. 

The predictions for Raa from the various models at 
y/SNN = 2.76 TeV turn out to be quite dramatically dif- 
ferent in normalization, shape and their expectation for 
even larger Px- In particular, the differences between 
models are significantly larger than the statistical errors 
of the measurement, thus allowing for a clean discrimina- 
tion if the systematic errors (such as those with the p-p 
baseline) can be understood. 

ASW and YaJEM in the default mode predict a rather 
slow rise of Raa with Pt along with a strong suppression. 
This is not in agreement with the data shown for the 
default p-p reference, but would agree better with an 
alternative NLO scaled p-p reference (see [l| for details). 

Both parametrized elastic scenarios reproduce the 
shape of Raa with the default reference well, they mainly 
differ in the normalization. The difference between the 



PbPb 2.76 ATeV, 0-5% centrality 



• ALICE data 




*0 10 20 30 40 50 60 

Pt [GeV] 



FIG. 3: (Color online) The nuclear suppression factor Raa 
in 0-5% central -^sjvjv = 2.76 TeV Pb-Pb collisions com- 
puted in various models for the parton-medium interaction 
(see text) with model parameters adjusted to describe 0-5% 
central ^/snn = 200 GeV Au-Au collisions compared with 
the ALICE data d . 

parametrized elastic and the ASW model can be under- 
stood as follows: At RHIC conditions, only a narrow 
region of P{AE) around zero is effectively probed due 
to the steeply falling spectrum, as even a small shift in 
parton energy is equivalent to a massive suppression. In 
the parametrized elastic models, most of the weight of 
{P{AE))taa is contained in the region between zero and 
- 30 GeV (see e.g. Fig. 1 right in (Hi), for small Pose 
even more weight is contained close to the origin. This 
is very different for {P{AE))taa c omp uted in the ASW 
model — here (see e.g. Fig. 5 in [40]) the distribution 
is much flatter and contains a sizeable weight out to 100 
GeV. 

To illustrate this in more detail: In a schematic model 
neglecting (among other things) hadronization, Raa can 
be understood from the ratio of modified over unmodified 
parton spectrum, where the modified parton spectrum at 
a given px is determined by the number of partons es- 
caping without energy loss plus the number of partons 
available in the spectrum at pt + AE times the probabil- 
ity P{AE) of a shift by AE. If we assume a power law 
p^" for the parton spectrum, 

Raa ~ J dAE{P{AE))TAA f 1 + — ) • (23) 

It is evident from the expression that Raa at a given 
Pt is equal to the transmission term of zero energy loss 
plus a contribution which is proportional to the integral 
of {P{AE))taa from zero up to the energy scale Smax of 
the parton (since a parton cannot lose more energy than 
it originally has), seen through the filter of the steeply 
falling spectrum. Raa then generically grows with pT 
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since -Emax grows linearly with px, and the speed of 
growth depends on the weight of {P{AE))taa the re- 
gion from zero to -Emax and on the power n of the parton 
spectrum. 

At the LHC, more of {P{/S.E))taa accessible com- 
pared with RHIC due to the harder underlying parton 
spectrum, i.e. the power n of the 'filter' by which large 
/\.E in {P{AE))taa are suppressed is reduced. This 
translates into a stronger rise of Raa with Pt as 
compared to RHIC kinematics, and this rise is most pro- 
nounced for models where there is substantial weight of 
the energy loss probability density close to the origin. 
Thus, the parametrized elastic model with small Pcsc 
shows the strongest rise, while ASW shows very weak 
Pt dependence. 

Finally, YaJEM-D predicts the strongest rise of Raa 
with Pt as a consequence of Eq. pT]) . In addition to 
the rise expected from the way the effective energy loss 
probability is probed as outlined above, YaJEM-D thus 
contains an explicit mechanism which introduces a strong 
energy dependence into the MMFF itself. While the nor- 
malization of the curve falls below the default baseline 
data, the shape is well reproduced. 



PbPb 2.76 ATeV, 0-5% centrality 
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FIG. 4: (Color online) The nuclear suppression factor Raa in 
0-5% central ^sjvjv ~ 2.76 TeV Pb-Pb collisions computed in 
various models for the parton-medium interaction (see text) 
with model parameters refit to ALICE data [Ij]. To indicate 
the magnitude of the uncertainty, the alternative p-p reference 
result which is not shown with errors by ALICE has been 
given a 10% error band. 



C. Refit to data 

As discussed in the context of Fig. [5] above, there 
can be residual uncertainties in extrapolating the QCD 
matter fluid dynamics description from RHIC to the 
LHC. Therefore, in the second step, we allow a refit of 
parton-medium interaction model parameters to the AL- 
ICE data. We introduce a parameter R to quantify the 
amount of refitting which is needed, where R stands for 
the ratio of modified over unmodified parameter -ftTmcd- 
Here we would consider, say, a 25% change in the model 
parameter reasonable within the uncertainties of the hy- 
drodynamical extrapolation. The results are shown in 
Fig.H 

Following this procedure, we find that the ASW model 
can be brought into a rough agreement with the alter- 
native p-p reference data, but even allowing for a sub- 
stantial parameter readjustment it does not agree with 
the default reference data — the Pt dependence is too 
weak. YaJEM likewise follows the trend of the alterna- 
tive p-p reference data well. YaJEM-D on the other hand 
can be brought into good agreement in both shape and 
normalization with the default p-p reference data. 

Unfortunately, even provided that we are willing to dis- 
card a model based on deviations R < 0.75 or i? > 1.25 in 
the refitting procedure, we can at present make no such 
statement from the ALICE data due to the uncertainty 
in the p-p baseline. This stresses the importance of hav- 
ing a measured baseline. However, looking at Fig. [31 we 
see reason to conclude that already a larger range in Pt 
will allow to decisively rule out some models based on 
the wrong shape of the Pt dependence. 



V. DISCUSSION 

In this work, we have simultaneously analyzed the low- 
Pt and high-P^ spectrum of charged hadrons in central 
Pb-Pb collisions at y'sjvTV = 2.76 TeV. For hydrodynam- 
ical modeling describing the Iow-Pt spectrum, we have 
computed the initial conditions based on pQCD mini- 
jet production and saturation [2^ (2^. To account for 
the NLO and higher order corrections in minijet produc- 
tion, we have made an iterative fit to the measured LHC 
charged-hadron multiplicity. The main outcome of this 
computation are the produced transverse energy and the 
early formation time, Ti ss 0.12 fm, of the pQCD minijet 
system. We assume that is also the starting time tq for 
the hydrodynamical evolution but in charting the uncer- 
tainties of our modeling, we have shown the sensitivity 
of the computed Pt spectrum to tq. Also the sensitivity 
to the decoupling temperature, multiplicity fitting, and 
transverse profile of the initial energy density is explicitly 
shown. 

In our hydrodynamic framework, the uncertainty re- 
lated to the initial transverse profile can be considered to 
be the main uncertainty in the extrapolation from RHIC 
to LHC. We have for simplicity considered two possibil- 
ities here, the eBC and the eWN profiles, thus essen- 
tially ignoring the QCD dynamics that may cause the 
profile to change from RHIC to LHC. The profile should, 
however, be computable in the pQCD-|-saturation ap- 
proach [27| but since we do not have the needed NLO 
pQCD elements fully at hand yet [6^1, we leave this as 
interesting future work. We nevertheless can see that 
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within the uncertainties charted, the pQCD + satura- 
tion + (ideal) hydrodynamics framework works quite weh 
in reproducing the charged- hadron Pt spectrum up to 
Pt ^ 4: GeV, and that we obtain the agreement essen- 
tially without further tuning of the hydrodynamic pa- 
rameters from RHIC to LHC. Also interestingly, we have 
observed (see Fig. [1]) that once parton energy losses have 
been accounted for, the computed pQCD tail of high-Pr 
hadron production starts to dominate over the hydrody- 
namic component at Pt ^ 5 GeV. Comparison with the 
LHC data shows that the matching of these two compo- 
nents is very efficient in that it leaves fairly little room for 
hadron production components in addition to the hydro- 
dynamics and pQCD+encrgy loss in the cross-over region 
= 4 - 5 GeV. 

After getting the obtained hydrodynamical evolution 
of the background QCD matter in control, we have pro- 
ceeded to analyze the high-Pp part of the charged- hadron 
Pt spectrum. We have applied several models of parton- 
medium interactions, tuned to the nuclear suppression 
factor Raa at ^/sNN = 200 GeV Au-Au collisions, to 
y^SNN — 2.76 TeV Pb-Pb collisions where the primary 
parton spectra are significantly harder. We found that 
in principle the large kinematic lever-arm at the LHC 
translates into a significant power to distinguish various 
models, even given the uncertainties in extrapolating the 
bulk medium model to larger energy. Of particular im- 
portance here is the rise of Raa with Pt observed at 
y/sNN = 2.76 TeV which is intimately connected with 
the slope of the pQCD parton spectrum and reflects the 
way the energy loss probability distribution is probed. 

If one had no systematic uncertainty in the data at this 
point, two of the models we tested (ASW and YaJEM) 
could be ruled out already by the ALICE data obtained 
with the default p-p baseline. On the other hand, these 
models agree fairly well with the alternative p-p refer- 
ence data. However, in addition to the angular averaged 
suppression factor in central collisions considered here, 
other observables need to be studied. If we require that 



a model should also account for the observed spread be- 
tween the in-plane and out-of-plane hard hadron emission 
as observed at ^/snn = 200 GeV, then YaJEM would be 
strongly disfavoured 0| and YaJEM-D would be more 
consistent with the data than ASW. 

The fact that the simple parametrized elastic energy 
loss model which is known to fail for pathlength depen- 
dent observables is able to describe the scaling in y/s from 
RHIC to LHC rather well should be a stern warning that 
^/s and Pt dependence only probe particular aspects of 
energy loss models, and that agreement with a subset 
of available observables may not be enough to judge the 
validity of a model. 

As shown in Fig. [3l extending the measurements of 
Raa out to larger values of Pt at the LHC will pro- 
vide strong constraints and viability tests for the parton- 
medium interaction models. A combined analysis of Raa 
in the dependence on Pt , ^/snn, impact parameter b and 
reaction plane angle will be highly discriminating be- 
tween the available models even without having to resort 
to other high Pt observables such as triggered correlation 
measurements or fully reconstructed jets. Constraining 
the nature of the parton-medium interaction by leading 
hadron production in this way is thus the first step to- 
wards tackling the more difficult task of understanding 
the complete dynamics of a parton shower in the medium. 
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